library(data.table)
library(ggplot2)
library(here)
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data.csv")
text.size = 10
# agriculture
#############################################################
library(data.table)
library(ggplot2)
library(here)
dir = here()
wd = gsub("figures", "data/world bank", dir)
setwd(wd)
d = fread("data.csv")
text.size = 10
# agriculture
#############################################################
agriculure = ggplot() +
#geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.1) +
#geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_path(data = d, aes(x = agriculture, y = energy, col = country_code)) +
scale_x_log10("Service Sector (% of Employment)", breaks = c(2, 5, 10, 20, 50, 100)) +
scale_y_log10("C02 Intensity of Energy (kg C02 / kg oil equiv.)",
breaks = c(0.1, 1, 10 , 100)   ) +
ggtitle("go") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times")) +
annotate("text", x = 5, y = 0.8, label = CO2.r2, parse = T,family="Times", size=3)
agriculure
agriculure = ggplot() +
#geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.1) +
#geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_path(data = d, aes(x = agriculture, y = energy, col = country_code)) +
scale_x_log10("Service Sector (% of Employment)", breaks = c(2, 5, 10, 20, 50, 100)) +
scale_y_log10("C02 Intensity of Energy (kg C02 / kg oil equiv.)",
breaks = c(0.1, 1, 10 , 100)   ) +
ggtitle("go") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
agriculure
agriculure = ggplot() +
#geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.1) +
#geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_path(data = d, aes(x = agriculture, y = energy, col = country_code)) +
scale_x_log10("Service Sector (% of Employment)", breaks = c(2, 5, 10, 20, 50, 100)) +
scale_y_log10("C02 Intensity of Energy (kg C02 / kg oil equiv.)",
breaks = c(0.1, 1, 10 , 100)   ) +
ggtitle("go") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
agriculture
source('~/Desktop/service/figures/sector_plot.R')
agriculture
source('~/Desktop/service/figures/sector_plot.R')
agriculture
dir = here()
setwd("~/Desktop/service/data/world bank")
dir = here()
library(data.table)
library(here)
library(magrittr)
dir = here()
wd = paste(dir, "/raw", sep = "")
setwd(wd)
meta = fread("WDISeries.csv")
codes = fread("codes.csv")
d = fread("WDIData.csv", header = F)
info = d[1,] %>% t() %>% as.vector()
series.codes = d[,4] %>% t() %>% as.vector()
ids = which(series.codes %in% codes$code)
y = as.numeric(info)
y = y[-c(1:4)]
keep = d[ids, ]
sub = keep[V4 == codes$code[2]]
n = ncol(sub)-4
country = rep(sub$V1, n)
country.code = rep(sub$V3, n)
year = rep(y, nrow(sub))
i = 1
sub = keep[V4 == codes$code[i]]
values = sub[, 5:ncol(sub)]
long = values %>% t() %>% as.vector() %>% as.numeric()
output = NULL
result = data.table(country, country.code, year, long)
View(result)
output = data.table(country, country.code, year)
output = cbind(output, long)
View(output)
for(i in 1:nrow(codes)){
sub = keep[V4 == codes$code[i]]
values = sub[, 5:ncol(sub)]
long = values %>% t() %>% as.vector() %>% as.numeric()
output = cbind(output, long)
}
s.code = codes$code[i]
sub = keep[V4 == s.code]
name = names(output)
output = data.table(country, country.code, year)
s.code = codes$code[i]
sub = keep[V4 == s.code]
values = sub[, 5:ncol(sub)]
long = values %>% t() %>% as.vector() %>% as.numeric()
output = cbind(output, long)
name = names(output)
name
name = gsub("long", "s.code", name)
name
output = data.table(country, country.code, year)
for(i in 1:nrow(codes)){
s.code = codes$code[i]
sub = keep[V4 == s.code]
values = sub[, 5:ncol(sub)]
long = values %>% t() %>% as.vector() %>% as.numeric()
output = cbind(output, long)
name = names(output)
name = gsub("long", s.code, name)
}
output = data.table(country, country.code, year)
for(i in 1:nrow(codes)){
s.code = codes$code[i]
sub = keep[V4 == s.code]
values = sub[, 5:ncol(sub)]
long = values %>% t() %>% as.vector() %>% as.numeric()
output = cbind(output, long)
name = names(output)
name = gsub("long", s.code, name)
names(output) = name
}
View(output)
setwd(dir)
dir
fwrite(output, "data_format.csv")
source('~/Desktop/service/data/world bank/format.R')
setwd(wd)
country.only = fread("country_only_codes.csv")
d = d[V3 %in% country.only$Code]
d = fread("WDIData.csv", header = F)
